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thermodynwlc limit. Me focus entirely on electronic energies, and accord 
Ingly take the protons to be fixed at positions \k[). If the electrons 
are assigned position and momentum operators 1?^) and (p^J, then the 
Hamilton Ian for the system is 


♦ H ♦ m m 
ep pp 


(1) 


where 





(?a) 


(2b) 




(2c) 


(Mere m Is an electron's mass and e the Magnitude of Its charge.) 

The mem charge density of electrons In this systw Is p • -eh/Q, and 
the electrostatic energy of a unlfona charge distribution with this density 

Z 

Is I / drdr* 2 . The traditional Interacting electron gas problem 

2 o |f-r ' | 

examines the InergJ of an electron system relative to this term, 1*e., 
the corresponding Hamiltonian Is taken to be 


iiai 

i * i L 


u |V s I 


I / drdr' 


Pi 




The ground state energy of M has been studied extensively: 4 * 7 when 

* 1V3 .1 

expressed In terms of the linear density parameter r f • (30/4*N) w a Q 

(a 0 Is the Sohr radius), It Is approximately 


HE® (r $ ) • «(— - - 0.115 ♦ 0.031 In r f )«y («) 

r s 2 r s 

where the contributions are the familiar kinetic, exchange and correlation 
energy terms, ;he latter being the Nozl Ires -Pines approximation (see 
below). Note also that in the large system limit, assumed here, may 
be written 


H « 

f 2m 20 

eg 

where 


v c 00 

* «we 2 /k 2 

end where 

♦ ♦ 


1k*r. 




i v c (k>; e (k) ; # (-o 

C#0 


(5) 

( 4 ) 

(7) 


Is the one-particle density operator for electrons* 

If we take n(r) as the total charge density of the protons placed In 
compensating uniform background of negative charge, l.e.. 


n(r) 


p ♦ e 


l 4(r-R-), 

t 1 


(•) 


(3) 


then (1) may be regrouped as 
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H • ♦ 
« 9 


♦ I /* drdr* 

20 .♦ ♦ . 


The tern (9b) can now be recognized as the Madelung energy, H E ^( r s *( R A D* 
for a static array of Ions in a uni font compensating background. As 
Indicated, Eg depends upon the density and upon the crystal structure of 
the proton lattice. However, the structural dependence is known to be 
weak, 1 and Eg may be accurately estimated using the spherical cell 
approximation In which each proton is surrounded by a sphere just large 
enough that the net charge within It vanishes (such a sphere has radius 
r & a 0 ), and In which the overlap of spheres is Ignored. This neutral sphere 
has electrostatic energy 


Vo 2 


• /* W «-r Z dr, (1- < 7 ^> 3 ) 


which accurately approximates the structural result 


E « * IS & *c ( *> V*S<4) 
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(For crystalline lattices p p (k) - N where (K t ) Is the set of 

reciprocal lattice vectors.) notice that the Madelung part of the 
Hamiltonian (l.e., 9b) may be rewritten exactly as 



irfiere C t Is the unit cell associated with site t. The spherical cell 
approximation therefore calculates the Madelung energy by neglecting the 
last term of (13) and by replacing C Q with a sphere of equal volime. 

Next we examine (9c), which is the electron proton interaction energy 
relative to the Interaction energy between protons and a uniform compensat- 
ing background. We compare two methods for treating this term: first the 
structural expansion technique, which was applied to a point-ion system 
(metallic hydrogen) by Hammerberg and Ashcroft.* Second, the MignervSeltz 
method irfilch Is developed further here. To obtain the structural expansion 
we first recast (9) Into the form 



irtilch Is valid In the thermodynamic limit, iy introducing a coupling 
constant into the last term we may use a theorem of Pauli* to write the 
ground state energy as 



t 


*°(V(* t » * *2 9 < r s > ♦ 

•; I » c (k)i p (-f) | <»|; t (»)|xXU (IS) 

iticrc |x> represents the exact many electron ground state for H at coupling 
strength a. tote that (9c) depends on structure through p p (-k). The 
structural expansion calculates <*|i # (k)|x> by developing It In orders 10 
of v c (k)p p (k). The linear tens In response theory Is 

p(X)v e <k) ( 16 ) 

•** r * I* the first order polarliiblllty of the Interacting electron 

gas. At this level of approximation 

^ * 2S ( l 0 <”> 

Thus, the structural expansion explicitly uses the result (4) of the Inter- 
acting electron gas problem, and Its results may be systematically Improved 
by applying higher orders of response theory. 1 

In contrast, the tflgner-Seltz 1 method Involves *proxinat1ons that are 
less well defined, and It cannot directly draw upon the standard electron 
gas result (4). However, It lends Itself more readily to direct physical 
approximation: for example. It calculates approximate one-electron wave 
functions >d»ich ore Inaccessible via the structural expansion. (The 
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numerical results of the two methods are. In fact, quit* similar.) In the 
Wlgner-Seltz method, one first Investigates approximate single-body wave 
functions 4j(r ) by making an Initial assumption (to be Improved In Section 
IV) that when an electron's coordinate Is found within the tfigner-Seitz 
cell of a given proton, the repulsion of the H-l remaining electrons Is 
exactly canceled by the attraction of the N-l remaining protons. In this 
situation an electron is Influenced only by Its nearest neighbor proton 
and «£{r) then satisfies the single-particle Schrddlnger equation 

^ . !?, , {(?) . <». 


where r may now be restricted to a single Mlgner -Seitz cell. The boundary 
condition on this equation Is most readily incorporated by invoking the 
spherical cell approximation, in which 11 the £ • fi wave function Is 
required to be spherically symmetric and satisfies 

(dp (r)/dr) - 0. (20) 

s 0 

(When t # 6, the appropriate boundary conditions are more complex,* but 
these conditions will not be needed In the following.) Thus the first step 
In the Wigner-Seltz calculation Is to determine the It • 0 "bottom of the 
band" energy by solving a radially Sjmmetric one-body SchrMinger equation 
within the spherical unit cell. The total ground state energy for the full 
Hamiltonian (1) Is then estimated using the techniques described In Section 
III. 
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III. TOT At MUST: SPfCftlCAL CELL *9901 INITIO* 

If m assessed the I ypr oilitt one-electron mvc functions. ^(r), 
ehlcb arc occupied In UN ground state, Utn we could evaluate the exact 
HMlItonlan (1) In the Slater determinant constructed from these wave func- 
tions. The result Inf Hartree-Focfc-llfcc approximation for the ground state 
energy Is just 


E HF ' i J Y (f,( * s’ 2 *! ; y v ) * f ( ? > (21*5 

i.s b 1 * |' * « t | * 

•-7.! .! //^At^t-IVVI'IWI* <»> 

21 k.s k‘.** QB Jfj - r 2 | * 1 

C .»* oo 1 |?,-rj| - k 1 k 2 k- 2 k‘ 1 


♦ I I* 

2 tin ♦ ♦ 

I*. • M 


uhere the f.s sums extend oxer the occupied levels (s represents the 
electron spin projection). Me assume the wove functions are normal lied by 

/ Br |^(r)| 2 * I. (22) 

Each Integral over the volmne In (21a-b) may be broken Into a sum over 
lattice sites m plus an Integral over the corresponding IHgner* Seitz cell 
C^. In the spherical cell approxlmetlon all of the "cross terms" 
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, £- 1/2 

£1 ♦ i_5 ijs « o 

d, 2 » 4 


uhlch liids tr the solution, regular at the origin. 


(26) 


* 0 (r) - Ae** /2 A<1 -t' ln ;I.y) (22) 

where N(a;b;y) Is t u — er *s confluent h j perqeonetrlc function. 11 * 1 * 
Application of the boundary condition (25) then loads to the slnple non- 
linear eigenvalue equation 

(| t. - 1 - C l/2 W1 - £J l/2 ;2;jr s > ♦ U ♦ ^ 1/2 )N(-^ 1/2 ;2;jr t ) • 0, 

2 (28) 

Tb* binding energy t b It easily found by solving (28) numerically: in 

U« regime of Interest rummer's function may be evaluated to desired 
accuracy by simply retaining more and more terms of Its pouer series 
representation. The norma Hied uave function for r $ ■ 1.5 Is plotted 
In figure 1. 

The wave function (27) has a cusp at the origin, as do all the hydro- 
gen atom uave functions. In fact, all these cusps are Identical In that 


(6(1" ♦)/«*r) rK> • - «i*. 


( 29 ) 


Tols cusp condition has recently been proven to apply »mder very peneral 
conditions: 1% It util eld us In later calculations when no exact solution 



1 ? 



which Is once apaln Whittaker’s equation. The P(r) dependent quantity In 
(3!) Is given by 
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This tera cm be evaluated using the power series representation for 
Runner *s function, tut convergence Is slow enough that this Is not 
practical for r $ greater than about 3.5. For larger values of r s the 
Integral Is readily evaluated nunerlcally. The regaining terns In the 
Hartree %m can also be evaluated, but only with considerable effort. This 
•ay be avoided, however, by noting that (35) Is the largest tern in the 
Hartree sun. This follows fro* the observation that the low energy k • 0 
wave function Is Itself auch larger near the center of the cell (irfiere 
potential energy Is low) than at Its edge, whereas the high energy wave 
functions are wore uni for* In amplitude. On the other hand, a lower bound 
on the sum Is necessarily given by the Hartree Integral for a flat 
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wave function, namely 


-“4a, 


/ drr 2 / drV* = 

0 4*(r s i 0 r 0 .4v(r s a 0 ) 3 


At r & * 1.5, the upper bound (36) Is only 5% larger than the lower 
bound (37). This fortunate numerical coincidence laplies that little error 
Is Introduced by taking the sue (23b) as the average of the two bounds (35) 
and (37). This Is precisely the approx laat Ion that we will enploy. 

The reaalning exchange ter* (23c) Is traditionally difficult to 
evaluate. Wlgner and Huntington 1 argue that it Is closely approximated by 
the uni fora Interacting electron gas result -(0.916/r s )lly, and we also 

accept this result. 

0 

Finally, t hf oust be corrected by adding a correlation energy. We 
will again use a uni for* electron gas result, namely the Itozltres -Pines* 
interpolation foraula (-0.115 ♦ 0.031 In r s )fty. Our final approximation 
for the ground state energy per electron Is thus 


0.115 ♦ 0.031 in r s ]Hy, 

(38) 


where appears In (36). This result Is plotted In figure 2: the zero 
pressure density corresponds to r s • 1.55, at idilch point the cohesive 
energy (relative to wel 1 -separated electrons and protons) Is -1.078 fty. 
The results are nuaerlcally quite close to those of Chakravarty et aK ,* 
which were obtained using a completely different aethod. 


IS 


lecMM Ux Naelltonlan H explicitly spin -Ind ependen t , the Po*t 
dl icesslon Is readily extended to consider a spin polarized yround state. 
In tills case the kinetic onergy term 2.21 It * It cones 3.51/r^, and tke 
exchange tem -O.BU/r^ becomes -l.lS4/r $ . Otherwise equation ( 38) Is 
ween an pad. These results are also «raphod In figure 1: paramagnetic and 
spin aliened en erg ies are seen to be equal at • 3. SO. 

IV. SCftECftlK IFfICTS 

In Section 111 wt found approximate single -body wave functions from a 
Schrddlnper equation which Ignored electron-electron repulsion. 
Subsequently, the total enerfy was obtained by evaluating the exact 
Haul 1 ton Ian in a state constructed from these wave functions. Better 
wave functions, and nore accurate energies, nay be obtained by utilizing a 
scr een ed potential. The SchrSdlnger equation then becones (hereafter ue 
use atonic units elth By • a^ • e*/2 • 1) 

(-»*-£♦ IN?)) ,.(?) . c. f. <?). (39) 

uhe*e U(r) Is sene nreregt potential arising fren electron -electron 
repulsion. Me new exanlne two such scree^i*^ potentials: first the 
potntlal of a n o g live m%t charge unlfomly distributed throw-out the 
dlgner-Seitz sphere, namely 

v f > • l-s i*t 7 - r *)- 

r. 


( 40 ) 



for the ftrtt case the k ■ 0 SchrOdtofe * efwetlon cm he set red eaactly 
using the P oner series nethod. For the m cm* the SchrMlaftr efwet Sm H 
enelyted varletlonelly. 

He bee In *1 *>t'*9 thet the tetal Henlltenien (1) cm be trtelelly 
rewritten m: 


M • H ♦ J l N?,- l f ) * U (? t - I,). («2) 

mere U(r) Is now chosen te renlsh for r ewtstde the Wlfner-Selti cell. Ijr 
eeeluetlnf this Heal I ton ten In e Sleter dot eminent co n pesod of solutions 
to (39). we f'.nd thet the Hertree-feck-Hke enerfy (23) I* ottered only by 
the edd It Ion of the (nefett.e) tem 

.1 l ( e? |*.<r)J* «r). («) 

This tun cen be estt noted (es the Hertree um wet), by eeereflnf the lower 
bound 


-4o{ * r*dr * 0 2 (r) U(r) 


(«) 
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m4 t ft* #ftr «r* 1*9 *Hoo#m1ot «*r# f ooctlaw* 1o (43)) 

- -St- £‘ rV d(r). («») 

r* 

tm& lost express to* Is o*»ot to 

- 12/{Sr t ). (eat fere Mct{r«M detcattel) 

. * — ill— - JL ♦ e** 1 *^—- ♦ _i_ )J. (Tkoedt-feral eotcatlet). 
r,J 2 k TF *Tf *fr 

% tflt* to* ter 11 or Hortroo mi, tot toss o# occoroc/ lovolvoi 1o toil 
OoonOfot QrocoOort H aot sorlovt. toU tkot «t r % • 1.5 to* Im fcoo Oi 
41f for Of ool/ X (for tot* potootfoH). 

WO will troot too oolfona Uckfroot potootlol first. Tfco Oo wo d ory 
oolot pros loo for p^r) Is tot* 

(. -L . 1L.1.L- Or,* - r*»J* p (r) . C 0 a # (r). !«*> 

dr r, 

f«V«r) r * °- (47) 

AUkae* (id) I* wewtljf ael • *»teadard* er 'aeaed* differ eatUl 
tealilw, tt It (Ml If Mired rite* tke aetiwd ef freklalrt. ** Oaljr oae 
Mlettea I* refuler it Ut ertfta. ead tt It 



t 0 - 0. t, • I. % z • -1. *, • (l/d)(2-« # ♦ Vr,) 

*i " * * ll 0 ' T * *«-* * r* **-«** 


The ter let far i(r) (MMffit *wl«Ul| ead aatferaljr far ell r. HdMf 
taterckeafe differential tee ead t wet tee. Hence teaftff eeadtttee 1*7) 


0*1) * 0. <s *' 

i*tcii leedt te aa etfcarelae caret tea far aad *W It eeetljr ealead 
aaaertcilljr. ■elite Out tke cat* caadtttaa (29) H aaectljr tel lifted kjr 
acre tract lea (id). deteHe tfce aw af • ae a C e e leMlt pataadtal (let 
faataeU • of ref create 14). 

Tke effaettre aett aejr ke d eter ataad free aarettaat (II) ead (K). 
dare mm *(r) * -Z/r ♦ U,(r). Ike rater eat te let tea af (») 1* 
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Hr) • I 0< *t ’ l * 9 i * ' VZ 

'.-«raw 12 • 


the Hartrii. MCkatf*, mtt torrtltXion twri<n trt tr«»t«d a* they 
«tr« for • port CoalMb potcotUl. To tfco grow* «*• »<*»*»• <*> 

n Mtt tltoroforo add oar attlMta for tbo *M (4J). Tba cttaltlag onergy 
ottlMtc 1l plotttd agalatt r % lo ftgr« J. Tba alalMa ooorgjr, C • 

-1.0S2 «y. foil* ot r % - l.*5. 

Tho walytlt of too TOom* -F eral potaattal My be vproacbcd •*«*« the 
«mo ooaer terle* art hod. Mead. ae flai thot (4*) I* raplacad by 


*o * °* *i * l * % ‘ J, *■ * £ o *«-2 ) /v 


The cat* coadUlaa it agala tat lifted exactly. Hoaerer. bccaatc s, 

M p«ad * «po« ill **» fracedtay coefflctaat* * 0 to « (> j. ae hare aot been 
able U ettabllsh that the terle* act rally coarerge*. taerlcal erldeace 
tadlcetet that It doe*, la fact, cearcrge. bat oaly very tloaly. for r, 
greater thaa aboat 1.X, coarer grace I* tloa caoagh that the eaact goaer 
terle* tolatloa I* of little balp la practical calcaletloat. Thete are 
phyttcally I vrtad ralaet of r % . to a dlffereat approach It clearly 

refulrod. 

The varia U omI «1M proves to 9c an accurate for finding 

^ oitf g^(r). * Suita* 1c trial nova fcnctlon |^(r) oast satisfy the 



sjanctr/ condition 

♦ t (r f ♦ A) - O t (r t - a), (59) 

In order to assure that the boundary condition (20) Is satisfied in such a 
way that the wove function is saooth d r ■ r . In addition. It Is e*- 
pected to behave like the Is hydropenlc wove function near the trifle, 

l.e.p 

♦ t (r) - A e* r for r « r % . (54) 

loth of these refulranents are satisfied by the choice 



Dilt i t (r) does not awtonatlcaily satisfy the cos# condition (29). Mi 
rjey therefore adopt either of two variational strotofles: 1) vary a chosen 

paraneters In * t (r) iditle Inposlnf me constraint (nomallxatlon), or 
11) Increase the variational son by one term, vary o*l paraneters, and 
Inpor* two constraints (nomallxatlon, and the cusp condition). In either 
case the variational problon Is straightforward and the results are 
actually fwlte slallar. This technlfue can easily bn checked hr oslnf It 
to recalculate to the Gewlonb potential Itself. Iflth M - 5 In the 

s«, the variational and truncated series results differed by only parts 
per oiniofl throughout the ranpe 0 < r % < S. The r % • 1.5 variational 
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Mm function cannot be distinguished Um figure 1 because of the width of 
the tracing pen. 

The Her tree Integral (3$) and the "screening correction" Integral (44) 
are straightforward!/ computed fro* the k ' 0 variational wave function. 

The onl / regaining obstacle Is the effective mass, which can again be found 
b/ Integrating (32) numerically. Alternative!/, the band width (and hence 
m*) may be calculated free an Isotropic two plane wave node! appropriate 
for a spherical trllleuln zone. 1 * If K Is the Magnitude of the shortest 
reciprocal lattice vector (so that K/2 Is the radius of the first trlllouln 

tone and i vr^*. 1 * (K/2)* ■ (2V) 1 ), then In this Model the energy £(k) 

relative to i(£ • 6) Is given b/ 

el.. • t » l ‘ l ; 

Where € Q (k) * fc* and i(k) 1* the Fourier transform of the potential. Fron 
(St) we find an approximate hand width, 

6. - - |*(K)J - 6(0) • 5 (K/2) 1 . (57) 

which defines the effective Mass m*. The results of the two Methods are 
rather similar (the two effective masses differ b/ 53 at r % - 1.5). In 
the following we will use the results of the first Method. 

The ground state energy of this last model Is again plotted In figure 
3. The mlnlmw falls at r $ • 1.61, where I® ■ - 1. 038 8/, a value which 
Is In close agreement with the results of Twa and Mahan.* 



V. CONCIUSIOIIS 

It Is worth emphasizing that the work reported here Is not a statisti- 
cal mechanical calculation. Me do not exaalne all possible states of N and 
find which Is the most probable. Instead we have proceeded from the 
assumption that the most probable state for the density range of Interest 
Is a metallic one in which the protons occupy the sites of a travels 
lattice. The spherical cell approximation obviously eliminates any refer- 
ence to the choice of travels lattice, and also eliminates the possibility 
that the ground state Is some other configuration, such as a lattice with a 
diatomic basis (which Is the experimental structure at zero temperature and 
atmospheric pressure) or a "liquid" In tdilch the proton disorder Is due not 
to thermal agitation, but to quantum f Iwr+sdSluns . 17 It Is thus not sur- 
prising that our curves show physically Irrelevant negative pressure for 
r % greater than about 1.45. This simply reflects the well-known fact that 
at these densities the ground state is not a travals lattice metal, but is 
either a molecular solid or a metal In voexi stance with such a solid. To 
complete the equation of state for hydrogen, then, we need to find 
l°(r s ) for a lattice with a diatomic basis. w and c onnect the two curves 
using a Gibbs tangent construction. Of course, there remains always a 
possibility that a third configuration, with still lower energy. Is the 
true ground state. Accordingly our calculations do not profess to find the 
zero temperature equation of state for hydrogen, hut rather the ground 
state energy of electrons In the presence of a putative static travals 
lattice of protons. (The major difference between this model and a plau- 
sible picture of metallic hydrogen Is the lack of. Ionic motion.) The 
results do, however, show that If metallic hydrogen exists In a stable 
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or artist*)* state, then It eaist have a density greater than about 0.60 
gu/cn* (uhich corre s pond s to r % ■ 1.65). 
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Figure 1 
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